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We present a timing and spectral analysis of the variable X-ray source 3A 1954+3 19. Our analysis 
is mainly based on an outburst serendipitously observed during INTEGRAL Key Program obser- 
vations of the Cygnus region in 2008 fall and on the Swift/B AT longterm light curve. Previous 
observations, though sparse, have identified the source to be one of only nine known symbiotic 
X-ray binaries, i.e., systems composed of an accreting neutron star orbiting in a highly inhomo- 
geneous medium around an M-giant companion. The spectrum of 3A 1954+319 above > 20keV 
can be best described by a broken power law model. The extremely long pulse period of ~5.3 
hours is clearly visible in the INTEGRAL flSGRl light curve and confirmed through an epoch 
folding period search. Furthermore, the light curve allows us to determine a very strong spin 
up of -2 X 10~^iyii during the outburst. This spin up is confirmed by the pulse period evolution 
calculated from Swift/BAT data. The Swift/BAT data also show a long spin-down trend prior to 
the 2008 outburst, which is confirmed in archival INTEGRAL flSGRl data. We discuss possible 
accretion models and geometries allowing for the transfer of such large amounts of angular mo- 
mentum and investigate the harder spectrum of this outburst compared to previously published 
results. 
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1. Introduction 

3A 1954+319 belongs to a small class of Low Mass X-ray Binaries (LMXB), in which a neu- 
tron star orbits around an M-type companion The orbits are typically wide, making accretion 
via common Roche-Lobe overflow impossible. Instead, material from the wind of the companion 
is accreted onto the neutron star and thereby converted into X-rays. This accretion mechanism is 
very common among High Mass X-ray Binaries and their O- and B-type companions with strong 
stellar winds, but not well investigated for late type binary systems. In the case of 3A 1954-1-319 
the accretion stream onto the neutron star seems to be highly variable as irregular X-ray flares are 
observed as well as times where the source is hardly detected at all. Binaries hosting a neutron 
star and an M-type companion are commonly referred to as Symbiotic X-ray Binaries (SyXB), in 
analogy of Symbiotic Stars consisting of a White Dwarf and a late-type star. These systems were 
named after their optical spectra in which both components, the hot White Dwarf and the cold com- 
panion are seen. Such features are not seen in SyXB, as the neutron star is invisible in the optical. 
A close look in X-rays is thus indispensable for a description of those systems. 

Clearly confirmed as SyXBs are only 7 sources today [^, two more are under dispute. Only 
very few of the SyXB, Uke GX 1-1-4 or 4U 1700-1-24 have already been studied in detail with dif- 
ferent X-ray instruments, including high-resolution CCD detectors, and much is still to be learned 
about these systems and the properties of the wind and the accretion mechanism. Being easily 
observable in X-rays SyXB provide a unique opportunity for studying the wind properties of late 
type stars, via X-ray variations of the emission lines and absoprtion column. 

3A 1954-H319 was detected in the late 1970s with Uhum and Ariel V [|3|, |l2|]. Q analyzed 



Ginga data and proposed that the system could be a HMXB. No more detailed studies were per- 
formed until [^] identified the companion as an M-star using an accurate Chandra -position, plac- 
ing the system in the group of SyXB. Soon thereafter, |[T]] discovered a ~5h period in Swift/BAT 
data. This lead to closer investigations by who presented the first broad-band spectrum of 
3A 1954-1-319, using INTEGRAL and BeppoSAX. More spectral analysis was presented by 0, 
who also showed that the spectrum is consistent with the typical spectrum expected from an ac- 
creting neutron star and can be modelled using an either a cutoff powerlaw model or a thermal 
Comptonization model. 

A strong decrease of the pulse period from 5.17 h to 5.09 h during an outburst in 2004/2005 
and possible physical explanations thereof was reported by [Q] in Swift/BAT data. The period itself 
is most likely the spin period of the neutron star, as it is much too short for the expected wide orbit 
of the system and no other significant periods have been found between 5 x 10"^ - 3.5 x 10"*^*^ s 
[§, It is also unlikely that the observed decrease in the period is due to the orbit of the binary 
system, as the mass function would then lead to impossible values of the mass of the companion, 
assuming a canonical mass of 1.4 M© for the neutron star. 

In 2008 November INTEGRAL serendipitously observed a strong outburst of 3A 1954-1-319, 
obtaining high quality lightcurves and spectra, which are presented here. 

2. Lightcurves and Timing 

The long-term daily lightcurve of 3A 1954-1-319 as measured with Swift/BAT is shown in 
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Figure 1: Left: Lightcurve as measured with Swift/BAT. The tickmarks above indicate times of INTEGRAL 
observations, with the blue tickmarks showing times when the source was significantly detected in INTE- 
GRAL/lSGm. Right: INTEGRAL /ISGRl 20-100 keV Ughtcurve in November 2008. The data were binned 
to Science Window resolution. The insets show zooms and a more detailed lightcurve with 100 s resolution. 
Overplotted in red is the average pulse profile. 



Fig. H, left. It is clearly seen that 3A 1954+319 shows irregular flaring with some flares being as 
bright as almost 100 mCrab and other episodes where it falls below the detection limit of Swift/BAT. 
Tickmarks above the lightcurve mark times when 3A 1954+319 was in the field of view of INTE- 
GRAL/lSGRl. A relatively dense sampling of INTEGRAL observations was performed during 
a flare in 2008 November when Key Programme observations of Cygnus X-1 (PI Wilms) were 
scheduled. As seen in the left panel of Fig. [I] this flare was one of the brightest since the moni- 
toring of Swift/BAT, with an average flux over the outburst of ~ 8.6 x 10"'°ergcm~^s"^ between 
20-100 keV, i.e., ~60mCrab and a peak flux of more than 180mCrab. We extracted INTEGRAL 
data of that outburst using the standard Off-line Science Analysis (OSA) 8.0 to obtain spectra 
and using ii_light, as distributed with OSA 7.0, to obtain lightcurves with a time resolution of 
100 s. The INTEGRAL flSGRl lightcurve of that period is shown in Fig. |l|, right. Even though 
3A 1954+319 was in an overall bright state during that time the INTEGRAL data show that the 
X-ray flux is still highly variable and varies by a factor of ~20. Using the epoch-folding technique 
[Q], a period of ~5.3h is clearly detected in the data. Compared to the results by [^, the pulse 
period is significantly longer, i.e., the neutron star has spun down between 2005 and 2008. 

As 10] found a strong spin up trend in the 2005 flare, we split the INTEGRAL lightcurve 
into three equally long parts and searched in each one for the pulse period separately. We found 
a clear spin up trend in these data, with the pulse period decreasing from 5.34 h in the beginning 
of the outburst to 5.26 h in the end. The results of the epoch-folding are shown in Fig. |^ From 
this analysis we determined the pulse ephemeris to be (during that outburst): T^^^ = 54782.6897 d 
(mean of observation), P = 5.3085 h, and P = -2x 10"^ ^h. To check for consistency, we folded the 
lightcurve with this ephemeris to obtain a pulse profile and compared it to the measured lightcurve. 
As seen in the insets of the right panel of Fig. [T] the pulse profile describes the lightcurve very well, 
in the beginning of the outburst as well as close to the end. This is also the first time that the pulsed 
flux of 3A 1954+319 is directly visible in a lightcurve, showing that 3A 1954+319 does not show 
large pulse-to-pulse variations. 
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Figure 2: Left: Results of the epoch folding analysis, with the lightcurve separated into three equally long 
parts. Additionally the overall result is shown in black and the range of the 2005 pulse period is marked in 
green. Right: Long term pulse period evolution, measured with Swift /BAT (black) and INTEGRAL jlSGKl 
(red). 

In order to put our INTEGRAL results into perspective, we updated the Swift/B AT period 
analysis of 3A 1954+319, again using the epoch-folding technique. Using all publicly available 
Swift/B AT lightcurve data, local period determinations were performed for 20 d long segments. 
Figure § shows all successful Swift/BAT period measurements, that is, all data in which a period 
could be significantly measured. As the Swift/BAT data a impaired by a low signal-to-noise ratio, 
we also analyzed more INTEGRAL /ISGRl archive data. Additionally to the three data points from 
the 2008 November outburst we could obtain 3 more measurements of the pulse period during 
times with extensive coverage of 3A 1954-1-319 with INTEGRAL. All data points, Swift/BAT and 
INTEGRAL /ISGRl, are plotted together in the right panel of Fig. ^ As can be clearly seen the 
measurements for both satellites agree very well. The pulse period is slowly decreasing during 
quiescence but increases strongly during prominent outbursts, as expected from standard accretion 
theory One such outburst was measured by [Q] in 2005, the other one was caught serendipi- 
tously by INTEGRAL in 2008 November. Even though in between these outbursts 3A 1954-1-319 
shows smaller flares, no measurable spin up is seen there. Correlating the luminosity with the 
change in pulse period led to no conclusive result, however, this may be due to the sparse sampling 
of the pulse period between outbursts. 

3. Spectral analysis 

The INTEGRAL observations present the unique possibility to gain knowledge about the ac- 
cretion physics during a strong flare of 3A 1954-1-319, which was accompanied by a strong spin up 
of the neutron star. Due to the fact that the source is close to INTEGRAL /ISGRVs detection limit 
during pulse minimum, we created a peak spectrum of only those 50% of data, which were mea- 
sured at phases close to the major peak of the pulse profile, thereby improving the signal-to-noise 
ratio. Comparing pulse profiles at different energies above 20 keV did not reveal any strong energy 
dependence with pulse phase, indicating that the spectrum is not sensitive to the selected pulse 
phase and that the obtained spectrum is very similar to an overall phase averaged spectrum. This 
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Figure 3: Top: Specttum taken with JEM-X and ISGRI, together with the best fit model. Bottom: Residuals 
of the best fit. 



Table 1: Best fit parameters of the bknpower model. 
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filtering resulted in the selection of 75 Science Windows (ScWs) with a pointing offset of < 10° to 
3A 1954+319, leading to a pulse peak spectrum with 445 ks exposure in ISGRI and 22 ks in JEM- 
X. ModeUing these JEM-X and ISGRI data between 3-100 keV with the pubhshed models lead 
to unacceptable descriptions only. A much better description in terms of ^^^^ (x^ = 44.14 with 43 
d.o.f.) was obtained, when using the bknpower model, a model consisting of two power-laws with 
different photon indices Fi and T2, with the change from one photon index to the other taking place 
at a specific break energy fbreak- It was not necessary to apply systematics for this fit. The model 
was attenuated by photo absorption at lower energies, expressed as equivalent column density of 
neutral hydrogen atoms A^h- Data and best fit model are shown in Fig. ^, and the parameters of the 
best fit model can be found in Tab. |l]. The required new description of the model compared to the 
older spectra leads to an overall harder spectrum, as no exponential cut-off is visible anymore. The 
soft energy continuum has a slope which is consistent with previous published results, see, e.g.. 



4. Outlook & Conclusions 



We presented a timing and spectral analysis of INTEGRAL data of the large outburst of 
3A I954-H319 in 2008 November. We confirm the unusually long spin period and measure a value 
of ~ 5.3 h. During the outburst we found that the pulse period is drastically increasing, with a spin 
up rate as high as P = -2 x lO*'* Yh- This rate is a magnitude higher than the one found by for 
the 2005 outburst. As these authors already stated, a change of this magnitude can not be due to the 
orbital motion, as it would lead to a mass-function for the companion of > 5 X IO^Mq. On the other 
hand, if we assume that the neutron star is spinning close to equilibrium, i.e., that the Alfven radius 
is equal to the Keplerian co-rotating radius, then a magnetic field of ~ lO'^ G would be required. 
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In the standard accretion picture of [Q], this assumption is justified as only close to equilibrium 
the torque can easily change direction. However, more sophisticated accretion models have been 



proposed for other sources, including retrograde spinning accretion disks [ |10| ] and spherical wind 
accretion. If no stable accretion disk forms and material and momentum is accreted from the stellar 
wind, a torque transfer of this magnitude might be possible (see Postnov et al. elsewhere in these 
proceedings). Continued monitoring of 3A 1954+319 is necessary in any case to clarify possible 
models and check if the source continues to spin down during quiet phases, while strongly spinning 
up during flares. 

The pulse peak spectrum of 3A 1954+319 we presented is distinctly harder than the previously 
published spectra, but is otherwise featureless. Especially, we find no indication for a Cyclotron 
Resonant Scattering Feature (CRSF), as seen in many other X-ray binaries, where the neutron star 
possesses a magnetic field > lO'^ G. Such a feature is the only direct way to measure the magnetic 
field strength in the X-ray producing region and could clarify the size of the accretion region. If 
3A 1954+319 really has a magnetic field on the order of 10^^ G, it is not suprising that we do not 
see any CRSF, as the fundamental line energy would be > 10 MeV, where the source is not visible 
anymore. 

For a more detailed presentation of the results see a forthcoming publication by Marcu et al., 
2011, in prep. 
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